Queueing process with excluded-volume effect 
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We introduce an extension of the M/M/l queueing process with a spatial structure and excluded- 
volume effect. The rule of particle hopping is the same as for the totally asymmetric simple exclusion 
process (TASEP). A stationary-state solution is constructed in a slightly arranged matrix product 
form of the open TASEP. We obtain the critical line that separates the parameter space depending 
on whether the model has the stationary state. We calculate the average length of the model and 
the number of particles and show the monotonicity of the probability of the length in the stationary 
state. We also consider a generalization of the model with backward hopping of particles allowed 
and an alternate joined system of the M/M/l queueing process and the open TASEP. 

PACS numbers: 02.50.-r, 05.70.Ln 



I. INTRODUCTION 

The queueing process is a typical example of a Markov 
process One of the simplest queueing processes is 
of the so-called M/M/l type, where the arrival of cus- 
tomers and their service obey the Poisson point process. 
This model's stationary state is the geometric distribu- 
tion that varies with the ratio of the arrival rate to the 
service rate. The M/M/l queueing process has no spatial 
structure and particles do not interact with each other. 

On the other hand, the asymmetric simple exclusion 
process (ASEP) on a one dimensional lattice is one of the 
simplest Markov processes with interacting particles Q. 
In the ASEP, each site can be occupied by at most one 
particle and each particle can hop to a nearest-neighbor 
site if it is empty. The ASEP admits exact analyses of 
non-equilibrium properties by the matrix product ansatz 
and the Bethe ansatz Q. The matrix product form of 
the stationary state was firstly found in the totally ASEP 
with open boundaries (open TASEP), where each particle 
enters at the left end, hops forward (rightward) in the 
bulk and exits at the right end Q. Similar results have 
been obtained in various generalized ASEPs and similar 
models in one dimension with both open and periodic 
boundary conditions [||. 

In this paper, we introduce an extension of the 
M/M/l queueing process on a semi-infinite chain with 
the excluded-volume effect (hard-core repulsion) as in 
the open TASEP. Each particle enters the chain at the 
left site next to the leftmost occupied site, hops and ex- 
its following the same rule as for the open TASEP. A 
stationary-state solution is given by a slightly arranged 
matrix product form of the open TASEP. The normaliza- 
tion constant is given by the generating function of that 
of the open TASEP. 

This paper is organized as follows. In Section 2, we 
briefly review the M/M/l queueing process and the open 
TASEP. In Section 3, we define the model. In Section 4, 
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FIG. 1: M/M/l queueing process. 



we find a stationary state of the model. This will be con- 
structed in a slightly modified matrix product form of the 
open TASEP. We obtain the critical line which separates 
the parameter space into the regions with and without 
the stationary state. The critical line will be written 
in terms of the stationary current of the open TASEP. 
We also calculate the average length of the system and 
the average number of particles on the assumption of the 
uniqueness of the stationary state. We also show the 
monotonicity of the probability of the length (i.e. the 
position of the leftmost particle). In Section 5, we gen- 
eralize the model by allowing particles to hop backward. 
It is fair to say that almost every calculation in Section 
4 and 5 will be performed by using known formulae in 
studies of the open TASEP and the open partially ASEP 
(PASEP). In Section 6, we introduce an alternate joined 
system of the queueing process and the open TASEP. 
Section 7 is devoted to the conclusion of this paper. 



II. REVIEW OF THE M/M/l QUEUEING 
PROCESS AND THE OPEN TASEP 



A. M/M/l queueing process 

Let us consider the simplest queueing process, i.e., the 
M/M/l queueing process, as in Fig. [T] where N denotes 
the number of particles. Particles enter the system with 
rate a and exit the system with rate (3. (Customers arrive 
at the queue with rate a and receive service with rate (3 
at one server.) The M/M/l queueing process does not 
have spatial structure and is characterized only by the 
number of particles. 

The system is encoded by a Markov process on the 
state space Z>o and governed by the following master 
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equation for the probability P(N) that the number of 
particles is N: 



dt 
d 

lit 



P(0) =/3P(l) - oP(0), 



(1) 



P(N) =aP(N - 1) + /3P(N + 1) - {a + P)P(N), 



(2) 



for JV £ N. The M/M/l queueing process is equivalent 
to a continuous-time random walk on Z> whose jump 
rates to right and left directions are a and (5, respectively, 
with reflection at 0. 

A unique stationary-state solution is easily obtained as 



TV 



N=0 



Z \(3 

In other words, the equation 



N 



(3) 



=/3P(l) - aP(0), (4) 
=aP(N - 1) + (3P{N + 1) - (a + (3)P(N) (N G N). 

(5) 

with X}jv=o = 1 imposed has the unique solution 

([3]). We should note, however, that the normalization 
constant Z does not always converge, in other words, 
the stationary state does not always exist. If a < (3, Z 
actually converges to /3/(/3 — a) and the system has the 
stationary state. (Otherwise, Z diverges and the system 
has no stationary state.) Thus the critical line is a — 
P. We can see the uniqueness of the stationary state 
recursively as follows: 



We can set P(0) = 



From (H]), we have P(l) 
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• The relation (0 implies that if we suppose 
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for JVeN, then 
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FIG. 2: TASEP with open boundaries. 

B. TASEP with open boundaries 

Let us consider an interacting particle system, the to- 
tally asymmetric simple exclusion process on the L-site 
chain with open boundaries (open TASEP), see Fig. [21 
Each site can be occupied by at most one particle. Each 
particle enters the chain at the left end with rate a if it is 
empty, hops to its right nearest-neighbor site in the bulk 
with rate p if it is empty, and exits at the right end with 
rate (3. Let us write ta = if the jth site is empty and 
Tj = 1 if it is occupied by a particle. The system is for- 
mulated by a Markov process on the state space {0, 1} L . 
The master equation on the probability P(ti, . . . ,tx) of 
finding a configuration (n, . . . , t£) is as follows: 



= a{2T 1 -l)P(0,T 2 ,...,T L ) 

E\ . _, 3 3+1 

[Tj+1 ~ Tj)P(Ti, . . . , 1, ,...,t l ) 

3=1 

+ /3(l-2r L )P(r 1 ,...r L _ 1 ,l). 



(10) 



In contrast to the M/M/l queueing process, the state 
space of the open TASEP is finite and thus it always 
has a stationary state. Moreover the open TASEP is 
irreducible, therefore the stationary state is unique. Der- 
rida et al found the stationary-state solution to the open 
TASEP in the following simple form 0]: 



P(n,...,T L ) 



Z L (a,[3,p) P 



X n \V(f)) 
(11) 



where Xo = E and X\ = D are matrices, (W(u)| and 
|V(w)) are row and column vectors, respectively, and 
Zi,(a, /3,p) is the normalization constant 



The average number of particles can be easily calcu- 
lated as 



(AOm/m/i = E NP ( N ) - J^Z- 

N=0 ' 



(8) 



In the stationary state, the current of particles through 
the server is nothing but the arrival rate a: 



(i J2 P(N) = a. 



(9) 



N=l 



Z L (a,l3,p) 



(W(f)\(D + E) L \V(£)). (12) 



The matrices and the vectors should satisfy the following 
relation so that the matrix product form (|11|) actually 
gives the stationary-state probability: 



DE = D + E, 

(W(u)\E = -(W(u)\, 
u 

\V{v))D = -\V{v)). 
v 



(13) 



Set (W(u)\V(v)) = 1 without loss of generality. The 
following representation satisfies the algebra (fH 
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(14) 



(15) 



where k = \/ {w + v — 1) /wv, a = (1 — w)/w and 6 = 
(1 — v)/v. Note that, in the original paper the bulk 
hopping rate p was set to be 1 and other representations 
of the matrices and vectors were found. For simplicity 
in what follows, however, we choose the representation 
(fT4| - ([T5| so that only the vectors depend on a, (3 and p. 
By using the algebraic relation (TT3"|h we can calculate the 
normalization constant as follows [4|: 



Z L (a,P,p) 



L 

£ 

0=0 



j(2L-j-l) 



\J+i 



(?) 



3+1 



L\(L-j)\ 



(16) 



The stationary current, for example, can be written in 
terms of the normalization constant Z^(a, f3,p) as 



J L (a,0,p) =p- 



(a,fi,p) 



Z L (a,P,p) 



(17) 



In the limit L — > oo, the phase diagram of the cur- 
rent J oc (a, (3,p) consists of three regions which are called 
the maximal-current (MC) phase, the low-density (LD) 
phase and the high-density (HD) phase (see Fig. [3|): 



Joo{a,j3,p) 



a/p) 
P/P) 



a,/3>£ (MC) 
a < min(/3, |) (LD) 
I (HD) 

(18) 



j3 < min(a, | 



The line a = (3 < | is called the coexistence line, where a 
shock between a low density segment and a high density 
segment exhibits a random walk Q: Joa{a, &,p) — a (l ~ 
a/p). 

Note that no particle can enter the chain if the leftmost 
site is occupied by another particle and thus the station- 
ary current is not equal to a. This means that the open 
TASEP is a "call-loss system." Recall that the station- 
ary current of the M/M/l queueing process, which is not 
a call- loss system, is a. 



III. MODEL 

Let us introduce a new model which is an extended 
M/M/l queueing process on a semi-infinite chain with 
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FIG. 3: Phase diagram of the open TASEP. 
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FIG. 4: Queueing process with the excluded-volume effect. 



the excluded-volume effect as in the TASEP. Figure 0] 
shows the model, where each site is numbered from right 
to left and L denotes the leftmost occupied site. Each site 
can be occupied by at most one particle. Each particle 
hops to its right nearest-neighbor site with rate p in the 
bulk if it is empty and exits with rate /3 at the right 
end. The rules of the bulk hopping and the exit are the 
same as for the open TASEP. However, the rule for the 
entering is different. Each particle enters the chain at the 
immediately left of the leftmost occupied site (or at site 
1 if there is no particle on the chain). One can put this 
model into the TASEP in the semi-infinite chain with a 
new boundary condition. 

The system is encoded by a Markov process on the 
state space 



5:=0U |J ({1}x{0,1} 



Let 



iL-l 



{0,(1), (1,0), (1,1), (1,0,0), 
(1,0,1), (1,1,0), (1,1,1), (1,0, 0,0),. 



(19) 



•}■ 



The leftmost particle of each element of S except can 
always be specified. Configurations such as (. . . , 1, 1, 1, 1) 
and (...,1,0, 1, 0, 1, 0) do not appear in S. 

Let P(l,T£_i, . . . ,ti) be the probability of finding a 
configuration (1, T£_i, . . . , T\) with the leftmost particle 
at .Lth site and P(0) be the probability of finding no 
particle on the chain. The master equation governing 
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the model is as follows: 

d „, 



dt 
d 

dt 



\ = /3P(1) — c\P(0), (20) 
P(l) = oP(0) +pP{l,0) - (a + /?)P(l), (21) 



= pP(l,0 ) r i _ 1 ,...,r 1 ) 

+ q;t l _iP(1, r L _ 2 , • . . ,n) - aP(l,T L _i, . . . ,n) 
+ p(T L _ 1 - l)P(l,0,n,_2...,ri) 

( t j ^ 7 j+i) "I 1 ! T L-ii • • • j i > o, . . . , nj 

3=1 

+ ^(l-2n)P(l,TL-i,...r 2 ,l). 



(22) 
For example, 

^P(l, 1,0, 1,0,1) 

= aP(l,0,l,0,l)+pP(l,0,l,0,l,0,l) (23) 
+ pP(l,l,l,0,0,l)+pP(l, 1,0, 1,1,0) 
- (a + 2p + /?)P(l,l,0,l,0,l), 

^P(l, 0,0, 1,1, 0,1,0) 

= pP(l,0,0,0,l,l,0,l,0)+pP(l, 0,1, 0,1, 0,1,0) 
+ pP(l,0,0,l,l,l,0,0) + /3P(l, 0,0, 1,1, 0,1,1) 
-(a + 3p)P(l,0,0,l,l,0,l,0). 

(24) 

We write the right-hand side of the master equation as 
(HP)(ti, . . . , tl) with the generator matrix H acting on 
the probability vector P = (P(0), P(l), P(l, 0), . . . ) T . 
Note that H is an infinite dimensional matrix. 



IV. STATIONARY STATE 

A. Matrix product form 

The problem is how to find the solution to HP = 0. 
We can see that a slightly arranged matrix product form 



P(0) 
P(l) 



1 



Z(a,P,p) } 

Vl \ A = 7( \ ^ (W(l)\V(f)), 



P(t l = l,...,n) 

Z{a,/3,p) 8p L 1 p 
(for L > 2) 



(25) 



gives a stationary-state solution. (This idea is applica- 
ble to a discrete-time version of the model 0.) Here 
Z(a,P 1 p) is the normalization constant which can be 
written as a special case of the generating function of 
the normalization constant of the open TASEP: 



oo r 

a L 



L=0 



oo r 



(26) 



with Z_i(p, P,p) = (3/p. See Appendix [Al for some sta- 
tionary probabilities calculated by using the algebraic re- 
lation (p~3|) . That the form (|25[) gives a stationary-state 
solution (i.e. HP = 0) can be proved by a similar cancel- 
ing to that for the open TASEP. Let us use a short-hand 
notation 



<w(i)| ■••!(!)) 



(27) 



Substituting the form (|23|) into HP(1, T£_i, . . . , ti) and 

o L — l 

multiplying it by pp L Z(a,/3,p), we find 



Bp 



L-l 



a 
a 



-Z(a,f3,p)(HP)(l,T L „ 1 ...,T 1 ) 



(28) 



=p-(EX TL _ 1 ---X Tl ) 
P 

P 

+ aT L - 1 -{X TL _ 2 ■ ■ X Tl ) - a{X TL _ 1 ■ ■■X TL ) 

a 

+ p{T L ^-\){EX TL _ 2 ---X Tl ) (29) 

L-2 . 

+ fa ~ r J+i)( X ^-i ••• D E ■■■X T1 ) 

3 = 1 

+ p(l-2 n )(X TL _ 1 ---X T2 D) 
=p{-(l-2T L ^){X TL _ 2 ---X Tl ) 

L-2 

+ ^(1 - 2r j+ i)(X TL _ 1 ■ ■■X T]+2 X Tj ■ ■ X Tl ) 

- - 2t j )(X Tl _ 1 ■ ■ ■ X T]+1 X Tj _ 1 ■ X Tl ) 

+ (1-2t 1 ){X tl _ 1 ---X T2 )} 
=0 (31) 

In going from (f2"5|) to ([3D]) , we applied the algebraic rela- 
tion CE3]): 

(...£>£;...)=(...£>. ..) + (.. .£...), (32) 
(E ■■■)=(■■ ■), p(...D)=p(...). (33) 

In our argument throughout this section, we assume 
that (j2"5)) is a unique stationary state of the model if 
Z(a,(3,p) converges and there is no stationary state if 
Z(a,(3,p) diverges. Recall that the stationary states of 
the usual M/M/l queueing process and the open TASEP 
are unique. 
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The following function will be useful in the next sec- 
tion: 







=l + J2^ L t(W(l)\(CD + E) 



L-l 



\V(§)) 



(34) 



L=l 



with a fugacity This can be regarded as a special case 
of the generating function of the normalization constant 
of the TASEP with a single defect particle (see section 4.3 
m @). The case where £ = a/p and ( = 1 corresponds 
to the normalization constant: 



Z(a,ft,p) = Z{pa/p,j3,p\ 1). 



(35) 




FIG. 5: Critical line 



Critical line 



fact, one can see 



The asymptotic behavior of Zl (p, ft-, p) as L — > oo is 
as follows [4|: 



Zl(p,P,p) 



4 
2-4 



^(^-) 0> E 

((i-Wp)Wp) ^ 



1-2/3/p 



(36) 



Thanks to this asymptotic form, we see that if the con- 
dition 



<*<! 



(/?>§) 
a</3(l-/3/p) (/3<f) 



(37) 



L=2 r 3 =0,l 



Z(a,ft,p) 

+ Z(a,P,p) f- 1 p 



I 00 L 

L=2 J 



(oo L—l ^ 



(40) 



C. Average values 



is satisfied, then Z(a,ft,p) converges. In other words, 
the critical line is given by 



a = a c 



E (ft > E ) 



(38) 



see Fig. O Note that the area (|37)) , where the normaliza- 
tion constant Z(a,ft,p) converges, is embedded in that 
of the usual M/M/l queucing process (a < ft). We re- 
mark that the critical line can be written in terms of the 
stationary current (fT5)) of the open TASEP with a = p: 



= Joo(p,/3,p). 



(39) 



Turning to the stationary current through the right 
end of the chain, we see that this must be the arrival 
rate a, because the model is not a call- loss system. In 



In this subsection and the next subsection, we assume 
that the condition (f37|) is satisfied. According to the 
formula (4.27) in [5j] we find 



= l 



where 



"=2 1 



i - (i + vera 
i - (i - vera 



In particular, 



Z(a,f3,p) = 



2ft 



2ft-p(l-r) 



(41) 

(42) 
(43) 

(44) 
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with r = a/1 — 4a /p. The average length (L) (the av- 
erage position of the leftmost particle) and the average 
number (N) of particles on the chain can be calculated 
by differentiating Z(p£, 0,p\ £) as 



(L) =e^lnZ(K,/3,p;l) 



2a /p 



(N) = —\nZ(a,0,p; C) 



t= a / P r(-l + r + 2/3/p)' 

(45) 

2a(l + r - (3 + r)a/p) 
C=i _ r(l + r)((l+r)0-2a)' 

(46) 



By the excluded- volume effect, these values are greater 
than the average length of the usual M/M/l queueing 
process: 

a 



(L) > (N) > (AOm/m/i 



— a 



(47) 



Actually, one can find 

(AO - (N) U /M/1 



2a 2 ((/? - a) + r(a + 0)) 



pr(l + r)(/3(l + r) - 2a) (0 - a) 



>0, 
(48) 

whereas the first inequality in (|47p is true by definition. 
(L) and (N) are expanded with respect to 1/p as 



(L) 



(N) 



a a 2 {20-a)l 
' 0-a + (0 -a) 2 p 



0((1/ P ) 2 



a 2 (3 1 



- + {(l/p) 2 ) 



(49) 
(50) 



— a [0 — a) 2 p 

and we can see a natural result that they approach 
{N)m/m/i m the usual-M/M/l-queuing-process limit 
p — > oo. Note that 

(L) = oo, (N) = oo (51) 

on the critical line a = ? and > | . 

D. Monotonicity of the length 

Let us consider the probability Al that the length is L 
(the leftmost particle is at site L): 



Al = 



1 



Z(a,0,p) 



TL-l, ■ 



l 



Z{a,0,p) \p 



(52) 



For L = 0, we set A = P(0). 

Thanks to the asymptotic form again (I25| . we can see 
that Al decays as L — > oo as 

' i(^BJ) 2 f^) L 0>R 

/^U \2f3-p) { p J ^ 2 



X, 



Z(a,0,p)0 



x < 



2 f 4a 



/?=§ • 



(53) 



When L is finite, Al possesses the property of the 
monotonicity with respect to L: 



< A 2 < Ai < A . 



(54) 



The rightmost inequality Ai < Ao is clearly true. We 
devote the rest of this subsection to the proof of Al+i < 
Al for L > 1. Let us use short-hand notations 



Z L ■= Z L (p,0,p), x:=^, 





(55) 



and the following alternate expression which can be ob- 
tained by transforming ()16D : 



Zl = 2_j CLL ' i x 

3=0 



,3 

3 * ' 



(56) 



where a LJ = jg^g^j ■ 



1. Case when < p/2 

Under the assumption a < a c = p/x(l — 1/a;), we find 
that 



1 / 1 



X \ X 



Cr 



-Z L - Z L -x < - 1 - - )Z L - Zl-i = t < 0, 



(57) 



where Cl = (L+r>!L! * s ^ e Catalan number. Thus, we 
have Al+i < Al- 



2. Case when > p/2 

The proof of the monotonicity for > p/2 (i.e. x < 2) 
will be somewhat more technical. Our goal is to show 
Z L < 4Z L -i, this implies that ^Z L < \Z L < Z L _ X and 

thus Al+i < Al- 

Before proving this in the general L > 1 case, we 
demonstrate it for L = 6: 

Z Cl = 132 + 132a; + 90a; 2 + 48a; 3 + 20a; 4 + 6a; 5 + x 6 

< 132 + 132a; + 90a; 2 + 48a; 3 + 20a; 4 + 8a; 5 

< 132 + 132a; + 90a; 2 + 48x 3 + 28a; 4 + 4a; 5 

< 132 + 132a; + 90a; 2 + 64a; 3 + 20a; 4 + 4a; 5 

< 132 + 132a; + 106a; 2 + 56a; 3 + 20a; 4 + 4a; 5 

< 168 + 168a; + 112a; 2 + 56a; 3 + 20a; 4 + 4a; 5 



= 4Z 5 



(58) 



where we used x 6 < 2x 5 , a; 5 < 2x 4 , etc. 

Let us go back to the general L > 1 case and intro- 
duce a sequence {&L,j}i<i<L-i defined by the following 
recursion relation: 



bh.j =2(&L,i+i - 4aL-ij+i) + ol. 



(59) 
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with aL—i,L — ^>L,L+i — 0. One can find that 
(2L-J-2)! 



(L+l)\(L-j)\ 

x (AjL 2 - 2(f - 4j - 3)L - j(j + 1)0' + 7))- 

(60) 



As long as bj^j > 4a L _ lj; 



(61) 



=b L ,j-ix j 1 +iaL-i,jX j . 

Let fc be an integer such that 

&L,fc < 4a L _i^ fc and 6 Lj > 4<z i _ 1 j (Vj > fc). (62) 

This is equivalent to 

k(k + 3) < 2L < (fc + l)(fc + 4) (63) 

and k is determined uniquely. For example, k = 2 for 
L — 6. Using the inequality (|6"Tj) repeatedly while j > k, 
we get 

•Zl = O£,o H 1- o-l,l-2X L ~ 2 + a L , L ^ix L ~ l + a LyL x L 

< a L ,a + ■■■ + a L ^ 2 x L - 2 + i>L,i-ia: i_1 

< ai,o H 1- b L . L - 2 x L - 2 + 4a L _i, L _i2; L ~ 1 

< • • • 

fc-1 L-l 



i=fe+i 



(64) 



The coefficients {aL.j}a<j<k-i in the first summation of 
the last line of (fM)) satisfy qlj < 4a,L-i,j'- 



< 



a L j - 4a L -i,j 

(j + l)(2L-j-2)l(j 2 +5j-6L) 

(L + l)\(L-j)\ 
(j + l)(2L-j-2)\(2j-AL) 



(65) 



(L + iy.(L-j)\ 



<o, 



where we used j 2 + 3j < k 2 + 3fc < 2L (see ([63]) ). 
Finally, we achieve 



fc-i 



L-l 



Z L < ^L-\,jX 2 + 4a L - hk x k + ^ ^L-l,jX j 



3=0 



j=k+i 



(66) 
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FIG. 6: Queueing process with exclusive hopping to both 
directions. 
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FIG. 7: PASEP with open boundaries. 

V. GENERALIZATION 

In this section, we generalize the model by allowing 
particles to hop backward with rate pq (q > 0), see Fig. 
O We will construct a stationary state and derive the 
critical line, arranging the matrix product form of the 
partially ASEP (PASEP) with the open boundary con- 
dition as in Fig. [7] 

The matrix product stationary state of the open 
PASEP is as follows: 

P(t u ...,t l ) 

=zj^: w ^--^\vo, (67) 

where the matrices Xq = E q and X\ = D q , the row 
vector (W(u;)| and the column vector |F(v)) satisfy 



D q E q -qE q D q = D q + E q , 

(W(w)\E q = -(W(w)\, 
w 

D q \V(v)) = -\V(v)), 



(68) 



and Zl(ck, f3,p,q) is the normalization constant: 

Z L (a,0,p,q) = ( W (f)\(D q + E q ) L \V(§)). (69) 

Some representations of the matrices and the vectors can 
be found in 0, H, The notation used in (f6"8"|) was 
chosen for convenience. However, we stress that, despite 
appearances, we have found no representation such that 
the matrices depend only on q and the vectors are inde- 
pendent of q. 

For q < 1, the normalization constant ZjJctJ3,p,q) 
can be written in the following integral form 



Z L (a,/3,p, q) 

(g,ab; q)^ f dz (z 2 , z~ 2 ; q)^ 
Aiti J K z (az,a/z,bz,b/z;q) c 



2 + z + z~ 



(70) 
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where a = -1,6 = - 1 and 

(xi,...,x n ' i q) oa = ni<„< m rio< i <oo( 1 - x n<f) is the 
g-shifted factorial. The contour K encloses poles z — 
a, qa, q 2 a, . . . and z — 6, qb, q 2 b, . . . , and excludes z — 
l/a,q/a,q 2 /a, . . . and z = l/b,q/b,q 2 /b, . . . . The 
asymptotic form of the normalizing constant of the open 
PASEP in the limit L — > oo has been obtained by ap- 
plying the saddle point method to the integral form ([70]) 
El: 



_ p(l-q) 



Z L {a,f3,p, q) 



' 4(ab;< ? ) 00 (g;g)^ 



2 ( 4 ) L 

2 f 4 \ L 

(a 2 ;g)oo / 2+a+a" 1 \ 

(b/a;?)^ V / 

( 



a, /3 > 



p(i-g) 



a 



= 



p(i-g) 

2 

p(l-g) 



</3 



< a 



&- 2 ;g)oo /^ 2+fc+fc- 1 \" 

/ 2+a+a -1 
V 1-9 



a = /3= 

a < min(/3, 
/3 < min(a, 



(a-l)(a :g)oo-L 
(a+l)^;?)^ 



a = /9 < 



p(i-g) 



(71) 



The stationary current of the open PASEP can be written 
in terms of the normalization constant as 



J L (a,j3,p,q) =p 



Z L ^(a,P,p, q) 



Z L (a,0,p, q) 
Noting the asymptotic form ([7T|) . we have 



(72) 



Joo(a,0,p,q) 



' p(i-g) 
a(l- 



(3(1 



p(i-g) 
g 

pa-?), 



a,/3> £^ 
/3 < m in(a, 

a < min(/3, silpl) 



(73) 



In the symmetric case q = 1, the normalization con- 
stant has the following simple form 1 101 ] : 

Z L (a, 0, p, q) = (7 + 1)(7 + 2) ■ ■ ■ (7 + L) (74) 

where 7 = ^ + ^ — 1. 

In the reverse-bias case q > 1, the normalization con- 
stant behaves in the limit L —> 00 as 



Z L (a,(3,p,q) ~ * 



(75) 



where A is a constant independent of L Q . 

Let us go back to the generalized queueing process with 
forward and backward hopping (Fig. [S]), which is gov- 



erned by the following master equation: 

|p(0)=/3P(l)-aP(0), (76) 
4p(1) = aP(0) + pP(l, 0) - (a + 13 + pg)P(l), (77) 

= pP(l,o,Tt_i, ...,n) -p<?P(i,Ti_ 1) ...,n) 
+ arL_iP(l,TL_ 2 , . .. ,n) - aP(l,rt-i, ■ ■ ■ ,n) 
+ p(rt_i - 1)P(1, 0, tl_ 2 • • • ,n) 

+ pt?(l - TL_l)P(l,TL_ 2 ...,Tl) 

( 7 "i- r i+i)"( 1 )TL-i,-.., 1 ,0,...,nj 

EJ+l J 
(Tj + 1 - Tj)P(l,T L _l, . . . , ,l,...,TlJ 

i=i 

+ /3(l-2r 1 )P(l ) r i _ 1 ,...r 2 ,l). 

(78) 

A stationary-state solution to this equation is given by 
the following form, which can be proved in the same way 
as (I58l-(I5I1): 



Z{a,f3,p,q) Z{a,(3,p,q) f3 

P(l,T L -x...,n) 

(79) 

where Z(a, /3,p,q) is the normalization constant: 

00 ^ 

Z(a,/?,p,g) = £ -^=I^-i{^,0,p,q) (80) 



L=0 



with Z-i(-£^,(3,p,q) = 13/ p. 

We now obtain the critical line assuming the unique- 
ness of the stationary state. In view of (|74|) and (|75jl . 
we find that the normalization constant Z(a, f3,p, q) con- 
verges only if the hopping ratio q < 1. Moreover, using 
the asymptotic form (|7Tj) . we find that the condition for 
the model to have the stationary state is 
for < q < |, 



a < a c 



a < a r = 



p(i-q) ^ > p(i-g) 



for ~ < q < 1, 



."( 



2(1-9) 
1 







P(l-?) 



/3< 



p(s+2-9) ' 



(81) 



(82) 



where s = y/ q(A — 3q). Note that the critical line a c is 
just the solution to 



(83) 



9 



(a) 







Q 



P 







a=a. c 






* 


p(l-?) 

o 
Z 




* 

✓ 

✓ 

✓ 

/ 










/* 

Jr 
A 


- 



p(l-flf) 



(b) 



/3 







a=a c 












✓ 

✓ 






/ 

/ 






/ 

/ 


p(s+2-g) 
2 


/ 


✓ 

/ 

/ 

/ 

■' 



p(s-q(3-2q)) 



a 



2(1-^) 

FIG. 8: Critical lines for (a) < g < | and (b) ± < g < 1. 
VI. ALTERNATE MODEL 



We introduce here an alternate joined system of the 
M/M/l queueing process and the open TASEP. This 
new system consists of a queue part and a TASEP part, 
see Fig. GO Each particle enters the system with rate a 
and joins the queue part. The queue part has no spatial 
structure, and is characterized by the number of parti- 
cles N. Each particle leaves the queue part and enters 
the TASEP part with rate a'. After entering the TASEP 
part, particles follow the same rule as in the usual open 
TASEP. This is a model of, for example, a production 
line with a material inventory. 

The state space of the Markov process encoding the 
model is Z>o x {0, 1} L . The master equation governing 
the probability P(N, n, . . . ,tl) of finding the configura- 
tion (JV,n, . . . , tl) is 



Q 



••• 











□ 







N 



L 



FIG. 9: Alternate joined system of the M/M/l queueing pro- 
cess and the open TASEP. 



Jt P(N,r u ...,r L ) 

= a(l-S m )P(N-l,n,...,T L ) 
-aP(N,n,...,T L ) 
+ a'r 1 P(A + l,0,r 2 ,...,l,0,...,r L ) 
-a'(l-n)P(N,0,T 2 ,...,l,0,...,T L ) 

L-l 

+ P^2( T 3 + 1 - Tj)P(N,Tl, ■ ■ .,Tj-tl,0,Tj+2, ...,Tt) 
3=1 

+ 0{1-2t l )P(N,t u ...t L -i,1). 

(84) 

For example, with L = 4, 



-P(0, 1, 1, 0, 1) =a'P(l, 0, 1, 0, 1) + pP(0, 1, 1, 1, 0) 



(a +p + /3)P(0, 1,1,0,1), 



(85) 



— P(5, 0, 1, 0, 0) =aP(4, 0, 1, 0, 0) -I- pP(5, 1, 0, 0, 0) 



+ /3P(5,0,1,0,1) 

- (a + a'+p)P(5,0,l,0,0). 



(86) 



It is difficult to find an exact stationary state of this 
model. One can expect, however, the critical line sepa- 
rating the parameter space to take the form 

a = J L (a',0,p), (87) 

because 

• In a stationary state, the current must be a. 

• If there is no stationary state and the queue part 
continues to grow, the TASEP part can be regarded 
as the open TASEP with a particle reservoir in the 
left end. Thus, the current must be Jl{oi' , f3,p) in 
this case. 

• These two values should be equal on the critical 
line. 

In fact, the critical lines for L = 1,2,3,4 with a' = p 
were calculated (although not rigorously) and found to 
agree with ([57]) [ll| . 



10 



VII. CONCLUSION 

We have studied an extension of the M/M/l queue- 
ing process on a semi-infinite chain with the excluded- 
volume effect as in the open TASEP. We found that a 
stationary-state solution is given by a slightly arranged 
matrix product form of the open TASEP and its normal- 
ization constant is given by the generating function of 
that of the open TASEP. The critical line which separates 
the parameter space into the regions with and without 
the stationary state is written in terms of the station- 
ary current of the open TASEP. We also calculated the 
average length of the system and the average number of 
particles. We also showed the monotonicity of the prob- 
ability of the length. These were derived by assuming 
the uniqueness of the stationary state. We generalized 



the model by allowing particles to hop backward and ob- 
tained its critical line. An alternate joined system of the 
queueing process and the open TASEP was introduced. 
We expect that its critical line can be also written in 
terms of the current of the open TASEP. 
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APPENDIX A: EXAMPLE 



where Z = Z(a, (3,p). These were calculated by using 
the algebraic relation ([To]) . 



The stationary probabilities (|25p for some configura- 
tions are listed here: 



ZP(l,0,l) = — , ZP(1,1,0) 



a 3 (p + /3) 



p/3 2 p 2 (3 2 
ZP(l,l,l) = p ZP(1,0,0,0) = ^, 

ZP(1, 0,0,1) 



ZP(1, 0,1,1) 
ZP(l, 1,0,1) 
ZP(l, 1,1,0) 



^ 2P(1,0,1,0) 


a%p + p) 
p 3 /3 2 


^ ^(1,1,0,0) = 


a 4 { P + 2(3) 
p 3 (3 2 


a 4 (p + /3) 




p 2 f3 3 ' 




a 4 (p 2 + P (3 + I3 2 ) 


ZP(1, 1,1,1) 



_1 

/3 4 ' 



